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Abstract 

The equations of state for neutron matter, strange and non-strange hadronic matter in a chiral SU(3) 
quark mean field model are applied in the study of slowly rotating neutron stars and hadronic stars. The 
radius, mass, moment of inertia, and other physical quantities are carefully examined. The effect of nucleon 
crust for the strange hadronic star is exhibited. Our results show the rotation can increase the maximum 
mass of compact stars significantly. For big enough mass of pulsar which can not be explained as strange 
hadronic star, the theoretical approaches to increase the maximum mass are addressed. 
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I. INTRODUCTION 



Stars are not only the focus of astronomers but also perfect laboratories for theoretical physicists, 
since extreme stars can supply conditions hard to be realized in earth-based laboratories. The 
neutron star was one of the most active actors on the stage of astronomy and physics since 1960s, 
when the first pulsar was observed [lj], and is still attracting more concentration. Many different 



aspects of neutron stars have been examined, such as the moment of inertia 



3|, the early stage 



of proto-neutron star jjj], the thermodynamic structure j^j, eh: However, because of the lack of 
the knowledge on the details of the matter at high density, the true equation of state (EoS) in such 
extreme object is still unclear. Great effort to determine the EoS of compact stars at extremely 



Mi 



or by 



high density have been made by many authors, either by the theoretical calculation 
observational constraint {s- 11]. At first, considering the extremely heavy mass and strong pressure, 
people regarded the pulsars as rotating neutron stars for a long time. However, after the pioneer 



work of Wit ten 



12| that the strange 



nucleus, the strange quark star 



13 



quark matter is comparatively more stable than the normal 



16J have attracted much attention (for a comprehensive review, 
see Ref. L7] ). I ~p to now, there is still no evidence strong enough to confirm the existence of strange 
quark star. The question is still open and need more careful study. 

In fact, the neutron star and the strange quark star are just very simple categories, and the 
real components in the star should be much more complicated. It is generally believed that the f3 
equilibrium can be achieved in the process of neutron-star formation and the hyperons will exist 
in neutron star. The consideration of hyperons in the EoS of compact star was back to as early as 
1960 



1 81 ] and got serious discussion in many later references 



compact 
0,0-21 



211 ] . The existence of hyperons 



can affect the EoS and consequently the mass-radius (M-R) relation remarkably. 



22 



231 ] . which 



In this paper, we will concentrate on a chiral SU(3) quark mean field model 
has been applied successfully to the study of nuclear matter, strange hadronic matter, nuclei and 
hypernuclei. In this model, quarks are confined within baryons and interact with mesons. The 
self-interaction of mesons is based on the SU(3) chiral symmetry. This model has been employed 



by Ref. 



24J to discuss the properties of static (non-)strange hadronic stars and neutron stars in (3 



equilibrium. However, the resulting maximum mass of the strange hadronic star is only lA5M sun , 
which is incapable to explain the existence of some massive stars, such as the recently confirmed 
data of PSR J1903+0327 with M = 1.67 ± 0MM sun However, we notice that the study 

in Ref. 24] limited to the static case. But in fact, rotation is a very general property of almost 



all stellar bodies, for example, PSR J1903+0327 has a very short rotating period 2.15ms. This 
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motivates us to extend the discussions of Ref. [24] to rotating cases and examine the rotation effect 



on the maximum mass. In this paper we will utilize Hartle's method [26] to study the slowly 



24J. We will also compare the 



rotating neutron star and hadronic star with the EoS given by Ref. 
M-R relation and other properties with those of the strange quark star obtained by the quark mass 
density- and temperature-dependent (QMDTD) model in Ref. [161]. 

The organization of this paper is as follows. In the next section, we will show the frame of 
Hartle's method applied in our calculation. In Sec. Ill, we will take the strange hadronic EoS of 
the chiral SU(3) quark mean field model to study the slowly rotating strange hadronic star. The 
effect of rotation on the maximum mass and moment of inertia will be discussed in detail. In 
Sec. IV, we will briefly show the numerical results of the non-strange hadronic and neutron star. 
Together with the results of the strange hadronic star in Sec. Ill and the strange quark star in 
Ref. [3], all these results for different stars will be compared. The last section is for summary and 
discussion. 

II. THEORETICAL FORMALISM FOR SLOWLY ROTATING STAR 

First we will present the formalism of Hartle's method on the slowly rotating star. This method 
treats the slow rotation as a small perturbation to the non-rotating structure. In static frame, the 
metric of a non-rotating, spherically symmetric star is given by 

ds 2 = -e u{r) dt 2 + e x{r) dr 2 + r 2 {d6 2 + sin 2 9d<p 2 ). (1) 

where e x ^ = (l — 2m ^ f m(r) = Jq Anr' 2 p(r')dr' . In hydrostatic equilibrium, the coefficients 
A and v are governed by the TOV equations: 

dv(r) 2 dP{r) 



dr p{r) + P{r) dr 

dP(r) [P(r) + p(r)] [m(r) + Airr 3 P(r)] 



dr r 2 



-y 2m(r) 



(2) 
(3) 



where P(r) and p(r) are the pressure and the energy density of the matter at radius r, respectively. 
The surface of the star r = R is defined by P(R) = and the stellar mass is M = m(R). Outside 
the star, the metric is simply of the Schwarzschild form, e u ^ = 1 — 2j ^-, so the inner and outer 
solution must connect at the surface of the star. These equations can be solved numerically by 
integrating outward from r = with a given center pressure P(0) = P c till the surface where P = 
is reached, then the main structure of a non-rotating star is obtained. 



Now we move ahead to the slowly rotating star with angular velocity 0(r, 9). The metric is 
not static and isotropic any longer, and the crossing term gte in the metric emerges. We will now 



261 ]. hereafter we only 



apply Hartle's method. Since this method had been shown in detail in Ref. 
write down the essential steps which are necessary in our numerical calculation. Following Hartle's 
expression, noticing Q is very small, the metric can be expanded up to the second order of 17: 

2m(r, 6)' 



ds 1 



e v{r) [l + 2h{r,9)]dt 2 + e 



A(r) 



1 + 



dr' 



r — 2m 

+r 2 [l + 2k{r, 9)}[d9 2 + sin 2 9(dip - u(r, 9)dtf] + 0(fl 3 ), 



(4) 



where h, m and k are corrections to the non-rotating metric of order Q, 2 , and u>(r, 9) is linear order 
of f2 in the gte component, which physically stands for the angular velocity of the local inertial 
frame. Define xu(r,9) = Q(r,9 ) — uj(r,9), the angular dependence can be studied by the vector 



spherical harmonics expansion 



27j. However, for the uniform rotating case f2(r, 9) = £1, one can 



find that w is a function of r alone as required by the boundary condition, so the equation obeyed 
by w is then simply 



26] 



1 d 

r 4 dr 



dw(r) 



dr 



+ -4 i 4) = o, 



r dr 



where j(r) = e 
j(r) = 1 and the solution is 



2 . The boundary condition at r = is 



dzu 
dr 



(5) 



0, while outside the star, 



(6) 



where J is the total angular momentum of the star. The moment of inertia is given by 

(7) 



J 

n' 



Besides the metric, the pressure is also affected by the rotation, which is denoted by an dimension- 
less parameter p*: 

p* = AP/(P + p). (8) 

The metric perturbation terms as well as p* can be expanded in spherical harmonics: 

h(r,9) = h (r) + h 2 P 2 (9) + ---, (9) 

m(r,9) = m (r) + m 2 P2(9) + ■ ■ ■ , (10) 

k(r,9) = k (r) + k 2 P 2 (9) + ■ ■ ■ , (11) 

P*(r,9)=p* (r)+P 2 (r)P 2 (9) + ---, (12) 



where the odd terms automatically vanish because of the symmetry. The field equation together 
with the TOV equations yields the equations for the I = correction terms: 

v? r " 5 r X ro . ( 13 ) 



drfiQ 
dr 


= 4irr 2 pQ(p -\ 


_ p) dp_ 
J dP 


dPo _ 


1 r 4 ?' 2 


fdw\ 2 


dr 


12 r — 2m 


{~d7j 



12 V dr J 3 dr 
+ H /r 3 jV\ 4vr(p + P)r 2 ^ m (l+8vrr 2 P) 
3 dr I r — 2m J r — 2m ^° (r — 2m) 2 ' 

where the boundary conditions are Pq(0) = m(0) = 0. Given a center rotating velocity w(0) = w c , 
the above formalism is enough to determine the changes on the stellar mass 5M and the mean 
radius SR, 

5M = m (R) + ^, (15) 
8R=- p * (R)p(R) / / d ^. (16) 

Higher order terms of the metric corrections are related to the configuration distortion away from 
an exact spheroid, which is neglected in our study. Under the linear approximation, we can use 
the relation w' c = w c j^ to obtain the desired value of w' c at the center from any initial value w Cj 
where £1' and f2 are the corresponding angular velocities. 

Finally we have to mention the so-called Kepler limit. With the increasing frequency, the 
tangent velocity of the matter on the stellar surface is also increasing. It is the gravitational force 
which keeps the matter from escaping. However, when the angular velocity is beyond a critical 
value, the gravitational attraction becomes insufficient to balance the centrifugal force, and the 
star becomes dynamically unstable. This onset is the Kepler limit given by the empirical formula 

VL K (M) w C(M/M sun )^ 2 (R/Wkm)- 3 / 2 , (17) 

where R corresponds to the radius of a non-rotating star with mass M, and C = 1.04kHz is a 
constant independent of the EoS [28|. When the rotating frequency approaches the Kepler limit, 
the stellar equator coincides with the innermost rotating orbit and the surface mass will shed away. 
Therefore we should always restrict our reasonable results within the Kepler limit. 



III. NUMERICAL RESULTS OF STRANGE HADRONIC STAR 

With Hartle's formalism, we can calculate the physical properties such as the mass, radius and 
angular momentum of a slowly rotating star with any given EoS. In this section we apply the 
EoS of strange hadronic matter from the chiral SU(3) quark mean field model [24| to examine the 
properties of slowly rotating strange hadronic star. 
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Before heading on to the numerical results, we should say a few words about the EoS at lower 
density. It is commonly believed that the neutron stars or quark stars are enveloped by a nucleon 
crust at the surface, where the density and pressure are much lower than in the center. The nucleon 
crust is supposed to be responsible for the glitch of the rotating frequency [29143 1|. and it can also 



affect the stellar radius significantly. To examine the property of the crust explicitly in the rotating 
hadronic star, we will compare the results of the strange hadronic star with and without the crust. 
The crust EoS from different models are generally similar to each other. In order to compare the 
rotating results with those of the non-rotating case given by Ref. [24], for ps < O.lfm -3 , we still 



use the same nuclear EoS as that of Ref. 



24j. namely, the EoS first given by Ref. 32], though many 
other good EoS for nuclear matter existed [y, [33J. At ps > O.lfm -3 , the EoS of strange hadronic 



matter shown in Fig. 4 of Ref. 



24| will be used. 



In FigHJwe show the M-R curves of crusted strange hadronic stars at different angular velocities. 
The curve turning over the peak value of M is kept with only a very short segment, since it is 
well known that the star becomes unstable at this segment. As can be clearly seen, when the 
rotation is imposed, the mass and radius are increased significantly. The faster the star rotates, 
the heavier and larger it can be. The round dot on each curve denotes the minimum stable mass at 
the corresponding angular velocity, which is obtained by the criterion of the Kepler limit following 
Eq. (I17p . Hereafter we will always use such round dots to denote the minimum stable mass on each 
curve, the segment with mass lower than the dot is unstable because of the violation on the Kepler 
limit. We also see from Fig|T]that the increase in the radius of the maximum mass at = 5000rad/s 
from the non-rotating case is about 3.55% only, so the shape deformation is indeed small. This 
confirms that our application of Hartle's method is self-consistent. In fact, as was demonstrated 
by Weber and Glendenning 34|, Hartle's method is applicable to stars with the rotation period up 



to 0.5ms, or in our language 0, = 12566rad/s, which is much larger than the fastest rotation speed 
observed yet. 

To illustrate the effect of the crust clearly, we also plot the corresponding results without crust 
in Fig J21 where the EoS of the chiral SU(3) model is still used when ps < O.lfm -3 . It is clearly 
shown that the shape of the M-R curve is dramatically changed. For a large range of M about 
0.1 ~ 1.2M sun , the M-R curves have different direction compared to the case with the crust in 
FigJTJ When the crust of lower density is considered, the radius of the star can be quite large when 
M is small, to which the outer crust contributes the most. 

Another important property about the crust is its moment of inertia, which is commonly used to 
explain the glitch of the rotating frequency. The glitch is a sudden change of the angular frequency, 
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FIG. 1: The M-R relation of strange hadronic star with crust at different angular velocities, where the 
round dots denote the minimum stable mass as required by the Kepler instability. 

which can be considered as the transfer of the angular momentum between the crust and inner 
part of the star. The recent observation requires that the crust must contain at least 1.4% of the 
total moment of inertia Now we calculate the ratio of the moment of inertia of the crust to 
that of the whole strange hadronic star. Since the moment of inertia can be calculated by 

I=^fl$V=9<W, (18) 
in our first order spherical approximation, it can be expressed as 

° I / 1 2m(r) Q, 1 ' 



The crust portion can be easily calculated by substituting the lower limit of the integration from 
to Rcrust-, which is determined by pB{Rcrust) = 0.1fm~ 3 . In Fig{3j we show the crust portion of the 
moment of inertia in the total star at different rotation speeds. We see that the crust contributes 
about 3.5% to the total moment of inertia at maximum mass, and this portion is larger at smaller 
mass. So the requirement of Ref.[9( is well satisfied. 

The most important thing we are interested in is the influence of the angular velocity on the 
maximum mass, which is shown in FigHJ where the curve of minimum stable radii vs. the angular 
velocity is also plotted. The minimum stable radius corresponding to the maximum mass can be 
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R/km 



FIG. 2: The same as FigQ] but the crust is not considered. 
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FIG. 3: The crust portion of the moment of inertia as a function of stellar mass at different rotation speeds. 

found directly from the M-R curve. We find that the minimum stable radius increases with Q. 
The dots denote the Kepler limit = 6454rad/s obtained from the formula in Eq. (|17p . It is the 
maximum rotation speed corresponding to the maximum mass of the non-rotating M-R curve. As 
can be seen from the M-R curves in Fig{H since the M-R curve of rotating star is higher than the 
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FIG. 4: The maximum mass and the minimum stable radius of strange hadronic star vs. the angular velocity. 

non-rotating case, the maximum mass on the M-R curve with such maximum rotation speed dose 
not lie on the boundary of Kepler limit, so this maximum mass can exist safely and lies far from 
the onset of dynamical instability. 

Besides the M-R curve, the I-M curve is also very important because it contains information 
about the inner property of compact stars. We show the I-M relation at different angular velocities 
in FigGD 



IV. COMPARISON BETWEEN DIFFERENT STARS 



In previous section we have shown the numerical results of the slowly rotating strange hadronic 
star in detail. Similarly, we can discuss the properties of the rotating non-strange hadronic star 
and neutron star. Here "non-strange" means that no hyperon is considered in the EoS, while 
"neutron" means only the neutron is considered. These EoS are also given in Fig. 4 of Ref. {24 1 . 
while the nucleon crust will always be considered from now on. The M-R curves for the non-strange 
hadronic star and the neutron star are shown in Figj6j and the I-M curves in FigJTJ respectively. 
All the Kepler limit points are dotted as in Sec. III. We see from these figures that the qualitative 
behavior is the same as the strange hadronic star. The masses and moments of inertia in these 
two cases are larger than the strange hadronic star, and the neutron star has the largest M and / 
among these three types. Another remarkable property shown in Figsj6] and [7] is that the curves 
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FIG. 5: The moment of inertia / vs. stellar mass M of strange hadronic star at different angular velocities. 




FIG. 6: The M-R curves of the non-strange hadronic star and the neutron star at different rotation speeds. 

of the same type are closer at bigger mass, while at smaller mass the curves of the same angular 
velocity are closer. This means that the rotational effect dominates in the lighter star with smaller 
central pressure. 

With all these results, now we can make some comparisons between different stars. They are 
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FIG. 7: The I-M curves of the non-strange hadronic star and the neutron star at different rotation speeds. 

the strange hadronic star, non-strange hadronic star and the neutron star obtained by the chiral 
SU(3) quark mean field model, as well as the strange quark star from the QMDTD model, whose 
properties are carefully studied in Ref. [r]] . 

A very important consistency check is that these stars have similar behavior according to the 
effect of the rotation. This is quite reasonable since the dynamical properties of a compact star 
are qualitatively similar. The centrifugal force from the rotation counteracting the gravitation can 
help the star bearing more mass and extending to larger radius. The details of these properties 
have been discussed clearly in Sec. III. 

One significant distinctness of the strange quark star is that its M-R curve has different direction 
from the other three types at small M, as was shown in Fig. 2 of Ref. [16]. The reason is that the 
crust is not considered in Ref. [3], so similar behavior is observed in our Figj2j When the star is 
very small, the crust thickness becomes important, which will change the shape of the M-R curve 
significantly. Generally, the nucleon crust can increase the radius by about 1km but has very little 
influence on the maximum mass. One difference between Fig. 2 of Ref. [16] and our Fig[2]is that the 
radius approaches to zero with decreasing M for the bare strange quark star, while it finally goes 
very large in our results of bare strange hadronic star. This difference can be understood when we 
pay attention to their EoS as the pressure approaching zero. In the QMDTD model, the energy 
density at zero pressure is finite, because of the confinement mechanism, so the pressure vanishes 
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FIG. 8: The curves of the maximum mass as functions of the angular velocity for the three types of stars. 

much faster than the strange hadronic matter. This property reduces the ability of such strange 
quark star to support a dilute outer layer as in strange hadronic star. 

Different EoS can quantitatively change the mass and radius of a star. The M-R curves of these 
stars show that, at the same condition, the strange quark star has the smallest radius, which is 
not changed by about 1km increase from the crust; while the neutron star has the largest radius 
and mass. The non-strange hadronic star has larger mass and radius than its strange relative. The 
maximum mass of strange quark star is a little smaller than that of the non-strange hadronic star. 
For example, at = 0, the neutron star has M max = 1.8M sun , for the non-strange hadronic star 
M max = 1.7M sun , and the strange hadronic star has M max = lA5M sun , while the strange quark 
star has approximately M max = 1.7M sun . These values are directly related to the EoS of different 
matter. In general, this property will not change: the steeper the P-p curve is, the larger mass 
the corresponding star will have. The rotation can not change the order of the physical quantities 
of these types of stars, either. To make this point clearer, we show the M max curves as functions 
of for different stars together in FigJHl Similarly, the R m in curves are shown together in FigJ^l 
where we see that the minimum stellar radius of neutron star is distinctly larger than the other 
two types. 
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FIG. 9: The curves of the minimum stable radius as functions of the angular velocity for the three types of 
stars. 

V. SUMMARY AND DISCUSSIONS 



In summary, employing the EoS of the chiral SU(3) quark mean field model and Hartle's method, 
we have studied the rotation effect on different properties of the strange hadronic star, non-strange 
hadronic star and neutron star. The most important effect of the rotation is that it can increase 
the maximum mass, which is listed clearly in Table HJ where O maiE is given by Eq.(17) according 
to the maximum mass of static star. We see that the range of the mass is significantly enlarged 
when the rotation is considered. However, the maximum mass of the strange hadronic star is 
still not large enough to explain the mass of PSR J1903+0327. As can be seen from FigJSJ the 
point representing the mass and the angular velocity of PSR J 1903+0327 lies well below the M max 
curves of the non-strange hadronic star and the neutron star but far above the curve of the strange 
hadronic star. So in the chiral SU(3) quark mean field model this star can not be explained as a 
strange hadronic star even the rotation effect is considered. 

We can easily understand the low mass of the strange hadronic star if we notice that its EoS 



35, 



3J, the 



is very soft. For other mean field hadronic model with hyperons, for example, Refs 
resulting maximum mass M max = 1.52M sun for the strange hadronic star and M max = \.&AM sun 
for the neutron star [35| at the static case are just a little larger than those we presented here. 
Their results still confirm PSR J1903+0327 is not a strange hadronic star, the same as our model. 
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TABLE I: The maximum mass of different stars with and without rotation 



type of star 


strange hadronic star 


non-strange hadronic star 


neutron star 


non-rotating M max 


lA5M sun 


1.7M sun 


1.8M sun 


M max at Qmax 


1.586M«,„ 


1.877M S „„ 


1.978M S „„ 


^max 


6454rad/s 


7039rad/s 


6968rad/s 



If we want to increase the M max of strange hadronic stars, the nonuniform rotation can be 
considered. As was pointed out in Ref . . the uniform rotation can just increase the maximum 
stellar mass by about 20% at most. Extend the discussion from uniform rotation to differential 
rotation, the upper limit of the stellar mass will increase remarkably 3814401]. Another possibility 



is to modify the soft EoS of strange hadronic matter, for example, the mixing phase of hadronic 
matter and quark matter. It has been proven that the EoS can be significantly steepened and 
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421 ] . In such case the 



the maximum mass be elevated when the mixing phase is considered 
deconfined quark becomes important and the interior part of the star might be made of quark 
matter, that is the so-called hybrid star 43j-|45|. It was shown that the hybrid star with a 2-flavor 



color-superconducting quark core can support very heavy star [461 ] . Another interesting question is 



471, 



48] . It is of interest to see the influence 



about the possibility of the "quarkyonic matter" phase 
of such phase on the M-R relation if they can be attained in the core of the compact star. However 
the discussion is beyond the chiral SU(3) mean field model we considered in this paper, because 
this model cannot provide the deconfinement mechanism. 
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